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Abstract This paper describes the development of a three-dimensional constitu¬ 
tive model for laminated veneer lumber (LVL) needed for new developments using 
this material. The LVL was manufactured in New Zealand from Radiata Pine. 
Experimental testing has been performed according to European timber testing 
standards. Block compression testing has resulted in modulus of elasticity values in 
the three material directions. Digital image correlation (DIC) technique has been 
used to determine the six Poisson’s ratios. Shear testing, whereby timber specimens 
were glued between two steel plates, has given stiffness values using DIC mea¬ 
surements. Experimental testing results have been compared with values found in 
literature. Results from this experimental testing programme have made it possible 
to create a three-dimensional elastic material model of LVL for the use in finite 
element analysis programmes. Although the material properties do not result in a 
symmetrical constitutive matrix, only minor adjustments are needed to gain the 
benefits of a symmetrical matrix. 


Introduction 

Timber is an orthotropic material, which means that material properties vary 
depending upon the orientation. Being a naturally grown product, there are 
numerous factors which influence the strength and stiffness, resulting in significant 
variability of material properties. Laminated veneer lumber (LVL) has the 
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advantage that defects are distributed, making the material properties relatively 
homogeneous compared to sawn timber (Buchanan 2007). 

New developments in the field of structural timber engineering, such as post- 
tensioned timber construction (Buchanan et al. 2011; Palermo et al. 2005) with 
complex connection behaviour (van Beerschoten et al. 2011a, b), require three- 
dimensional (3D) finite element modelling (FEM). To get reliable results from these 
models, it is necessary to have an accurate material definition to form the 
constitutive matrix. This material model consists of three moduli of elasticity ( E ), 
three Poisson’s ratios (v) and three shear moduli (G). 

The first publications on elastic properties in timber date back to the early 
twentieth century, of which an overview table is published in Hearmon (1948). 
Other publications on elastic constants of wood from similar time refer to these 
initial tests (U.S. Department of Agriculture Forest Products Laboratory 1955; 
Wangaard 1950). Later publications introduce different test methods for determi¬ 
nation of elastic properties. Gunnerson et al. (1973) described a plate testing method 
using two-way bending, which is concluded to be a good technique to determine 
elastic properties, but problems arise in calculation of Poisson’s ratios as small 
deflections cannot be measured accurately enough. Bodig and Goodman (1973) 
used this technique and compression testing. Their publication contains an extensive 
list of elastic properties of different timber species, including numerous pine 
species. A more recent publication (Bucur 2006) describes ultrasonic techniques 
which can be used in nondestructive testing. This method is compared with static 
compression testing by Goncalves et al. (201 1). The only values for LVL have been 
found in a publication by Janowiak et al. (2001), where several elastic properties of 
three types of LVL have been evaluated. An overview of elastic properties found in 
literature is provided in Table 1. 

Worldwide, there are a number of testing standards to determine mechanical 
properties of timber. These standards usually specify bending, compression or 
tensile tests for evaluation of modulus of elasticity. Shear modulus is often specified 
based on bending tests (single span or variable span), torsion tests and more recently 
shear field tests (EN 408 2010). No timber testing standard specifies tests for 
Poisson’s ratio, but ASTM El32-04 (2010) specifies a general test method for 
determination of Poisson’s ratio using extensometers during tension tests of 
structural materials. 

In literature, several methods for determining Poisson’s ratios can be found. Early 
publications on Poisson’s ratios in timber (Hearmon 1948) do not describe exact 
testing procedures except that testing was carried out using large-scale testing 
machinery. Sliker (1972) described a method using small strips of timber (32 
in. x 3.5 in. x 0.25 in.) in a testing machine with bonded electrical-resistance strain 
gages mounted in a rectangular pattern in the centre of the specimen. Results were 
consistent with a small variation, but only two material directions have been tested. 
Zink et al. (1997) published the use of digital image correlation (DIC) technique using 
white light speckle technique on compression specimens. This technique allowed 
multiple measurements spanning the entire specimen rather than from one or two 
points per specimen. This technique proved successful and the authors discovered that 
the Poisson’s ratios are not constant during testing, but decrease with increasing load. 



Table 1 Timber material properties from literature and experimental testing 
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A similar technique was used by Ling et al. (2009), whereby stochastic neural 
networks were used to approximate the displacement profiles. Niemz and Caduff 
(2008) have determined the Poisson’s ratios of spruce using a specially fabricated 
testing machine, capable of measuring transversal deformations with an accuracy of 
0.0015 mm. They tested 20 samples for each of the three different orientations. Results 
were in line with literature values, and variations between 17 and 62 % were found 
between measurements. Another publication (Garab et al. 2010) used DIC for 
determining the Poisson’s ratios under seven different growth ring angles, whereby a 
total of 182 specimens were tested. A video system recorded a 10 x 10 mm" section 
with a resolution of 950 x 950 pixels and used this to create full-field displacement 
profiles, resulting in low coefficient of variations (7-19 %). 

There have been numerous publications on shear tests focussing on shear strength 
of timber, i.e. Denzler and Glos (2007). The number of papers focussing on shear 
stiffness is limited, though several testing methods are published. Vibration time 
measurements have been used to determine the shear modulus of glulam beams 
(Gorlacher and Kiirth 1994). Divos et al. (1998) compared a static three-point 
deflection method with vibrational methods and concluded that static and torsional 
vibration methods are good, but the vibration method is more precise. A new block 
shear test (Sretenovic et al. 2004), whereby the specimen was glued between two 
sections of beech, gave promising results and, in comparison with ASTM D143 
2009 block test, showed a more uniform shear field. Another method, similar to the 
ASTM test setup, was proposed by Ukyo et al. (2010), and the DIC technique has 
been used to evaluate the shear modulus. A similar technique (digital speckle 
photography) has been used by Hassel et al. (2009), but this time on small square 
blocks, making it possible to test different material orientations, although only the 
rolling shear modulus has been tested. Brandner et al. (2008) used shear field 
measurements during standard four-point bending tests. This procedure integrates 
determination of shear modulus with standard testing methods. For this test method, 
large specimen sizes are needed, which is not a problem when glulam is used, but is 
problematic when analysing shear moduli in different material orientations. Dahl 
and Malo (2009) tested shear properties of Norway spruce using the Arcan shear test 
of small specimens in six different directions. Video extensometry was used to track 
displacements of a grid of dots on each specimen. Their test setup meant that they 
had to perform FEM analysis in order to derive modification factors for the shear 
modulus as shear stresses were not constant in the measurement area. The method 
published in this paper combines a standard testing method using small specimens 
with the use of the DIC technique to provide an easy way to measure the shear 
modulus in different material orientations. FEM analysis has resulted in a single 
factor to take the nonlinear stress distribution into account. 

Digital image correlation (DIC) is a technique for measuring deformations and 
strains (Pan et al. 2009). Software allows users to track displacements of points 
through a series of images taken during the experimental testing. This technique 
provides full-field displacements and strains, whereas traditionally used potenti¬ 
ometers and strain gauges only give point wise measurements. Although DIC is 
widely used in the field of experimental mechanics, there is still a great potential for 
this technique in the field of timber testing. 
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This paper presents experimental results of compression and shear testing on 
Radiata Pine LVL, performed according to European standards EN 14374 (2004) 
(requirements for LVL) and EN 408 (2010) (testing methods). The latter specifies 
pure block compression tests and shear tests with a relatively uniform shear 
distribution. DIC techniques have been used to determine the material properties 
needed for the formation of a 3D elastic material model. A full description of the 
compression strength and stiffness of LVL has been published in van Beerschoten 
et al. (2013), and this article provides a summary of these test results. 


Theory 


For the development of a 3D linear finite element material model, a total of nine 
material properties are needed for the constitutive equations; three moduli of 
elasticity (E), three Poisson’s ratios (v) and three shear moduli (G). These equations 
are shown in the form of the compliance matrix in Eq. 1. For LVL, having veneers 
rotary peeled and glued back together, there is no influence of growth rings and a 
clearly defined rectangular coordinate system can be used as shown in Fig. 1. The 
three material directions are (1) Longitudinal—parallel to grain; (2) Tangential— 
perpendicular to grain and parallel to glue lines; and (3) Radial—perpendicular to 
grain and perpendicular to glue lines. 
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For an orthotropic material like timber, there are a total of six Poisson’s ratios: 
v i 2 > v i 3 > v 2 u v 23 i V 3 U v 32 - Although it is commonly assumed that only three of these 
are independent, as finite element analysis programmes assume a symmetrical 
constitutive matrix. This means that off-diagonal terms have to be equal, which 
leads to relationships in Eq. 2. Validity of this assumption is evaluated in the 
discussion section of this paper. 


Vl 2 _ V 21 Vi 3 _ V 31 V 23 _ V 32 

E\ Ei ’ E\ £3 ’ E 2 E 23 



Materials and methods 

The LVL used for this research is manufactured from New Zealand Radiata Pine by 
Nelson Pine Industries Limited (2010). Specimens had a thickness of 45 mm and 
consisted of 13 veneers of approximately 3.5 mm thick. The material is specified as 
LVL11, with characteristic properties as shown in Table 2. These values are based 
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Fig. 1 Coordinate system and 
stress directions in LVL 



Table 2 Characteristic properties of LVL 11 (Nelson Pine Industries Limited 2010) 


Property 


Magn. 

Unit 

Modulus of elasticity 

E 

11 

GPa 

Bending strength 

/b 

48 

MPa 

Compr. strength par. to grain 

fc 

45 

MPa 

Compr. strength perp. to grain 

/c,90 

12 

MPa 

Tensile strength par. to grain 

/. 

30 

MPa 

Shear strength (edgewise, Gtl) 

/. 

6 

MPa 


on AS/NZS 4063:1992 (Standard New Zealand 1992) strength values are derived 
from the 5th percentile, corrected by an equation including the variability and 
sample size. Stiffness values are based on corrected mean values. Although values 
for one modulus of elasticity are given by manufacturers, this is based on four-point 
bending tests and include shear deformation (Nelson Pine Industries Limited 2010). 
Therefore, this value cannot be used for the constitutive equations. The timber was 
stored at ambient laboratory conditions, and density and moisture content were 
evaluated immediately following testing. Density was measured at an average of 

o 

581 kg/nr (CoV = 0.7 %) and moisture content at an average of 9.3 % (CoV = 
14 %) for shear testing and 8.9 % (CoV = 7 %) for compression testing. For shear 
testing, a two-part epoxy (Epi-Glue) was used to fix the steel plates to the timber 
specimens. Pressure was applied for about 4 h while the epoxy hardened after which 
specimens were placed in an oven at 35 °C for the epoxy to fully cure. 

Compression testing was performed using five or six replicates in the three 
different material directions for 45 mm thick LVL. European Standards have been 
used for compression testing. EN 14374 (2004) specifies requirements for LVL, 
although for testing it refers to methods outlined in EN 408 (2010). For parallel to 
grain stiffness, a small column with a height of ‘6 x dC (d = depth of specimen in 
mm) is specified. For compression perpendicular to grain, a block compression test 
is specified with a specimen height of 90 mm. Further specimen dimensions are 
given in Table 3. Testing was performed using an Instron testing frame with an in¬ 
line 150 kN load cell. Two small linear displacement potentiometers (10 mm travel) 
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Table 3 Specimen description and average dimensions for tested compression and shear specimens 


Test 

Parameters 

No. of tests 

Length (mm) 

Depth (mm) 

Height (mm) 

Compr. 1 

E 1 , V 12 , v1 3 

6 

44.7 

45.1 

271.2 

Compr. 2 

E 2 , V 21 , V 23 

5 

71.2 

45.4 

89.2 

Compr. 3 

£ 3 , V 31 ,V 32 

5 

71.6 

44.9 

91.1 

Shear 1 

G\2 

6 

295.0 

45.3 

44.9 

Shear 2 

G\3 

6 

272.7 

29.9 

44.8 

Shear 3 

G 21 

6 

296.8 

45.4 

45.3 

Shear 4 

G 23 

6 

272.9 

31.9 

42.6 

Shear 5 

G 31 

6 

295.3 

31.7 

43.9 

Shear 6 

G 32 

6 

295.7 

31.8 

45.2 


were placed on both sides at 20 and 80 % of the specimen height. Testing was 
stopped at 10 % strain deformation, as further deformation was deemed unrealistic. 
Stiffness is defined as the slope of load-deformation curve between 10 and 40 % of 
maximum load. 

Shear testing was also performed according to EN 408 (2010) with six specimens 
for each shear direction. Testing was performed in an Avery test frame with a 
capacity of 100 kN whereby specimens were loaded until failure. Specimen 
dimensions are given in Table 3. Due to material dimensions, specimens had a depth 
of 45 mm instead of 55 mm as specified in EN 408 (2010). Therefore, steel side 
plates were manufactured with a thickness of 15 mm, 5 mm thicker than specified 
by the standard, in order to keep the 14° angle of the specimen. Further details of the 
shear testing can be found in Dunbar et al. (2011). 

During compression and shear testing 18MegaPixel (5,184 x 3,456 pixels) 
RAW images of the specimen have been taken at intervals of 5 s using Digital 
SLR cameras, as shown in Fig. 2. These images have been analysed using 
software, a MATLAB script (developed by C. Eberl from Karlsruhe Institute of 
Technology, freely available for download from MATLAB Central), which was 
used to generate a grid of markers which were tracked through the series of 
images. For compression testing, 40-70 markers were placed between 20 and 
80 % of the specimen height. Vertical and horizontal displacements were 
converted to strains and averaged for all markers. For shear testing, grids of up 
to 280 markers have been used. 


Experimental testing 

Modulus of elasticity 

The modulus of elasticity ( E ) has been evaluated using the average displace¬ 
ment readings from potentiometers on both sides of the test specimen and load 
values from the load cell. Displacement values were divided by the distance 
between measuring points, and load values were divided by specimen cross- 
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sectional area. The modulus of elasticity has been evaluated using a trendline 
between 10 and 40 % of the maximum stress in the stress-strain graph. Further 
information about the compressive strength and stiffness can be found in van 
Beerschoten et al. (2013). 

Poisson’s ratio 

The Poisson’s ratios (Fig. 3a) of LVL have been evaluated from strain fields 
generated using DIC software. For each test, two series of images, one for each side, 
were taken using DSLR cameras. A virtual grid of markers was generated, as is 
shown in Fig. 3b. Each marker represents a square of 50 x 50 pixels which are 
tracked through the series of images. The coordinate of each marker is stored, and 
displacement relative to the first image is calculated. 

Average vertical and horizontal strains, calculated from displacements, 
during testing are shown in Fig. 3c. Also shown in the same figure is the 
Poisson’s ratio, with almost perfect linear behaviour after the first 5 pictures. 
The Poisson’s ratio has been based on an average value between image number 
30 and 60, whereby negative values were omitted. This range was different 
from the 10-40 % of maximum load which was used for evaluation of modulus 
of elasticity and shear modulus. During the first stage of testing, strains were 
too small to give reliable results and close to failure, the material behaviour was 
no longer linear elastic. For testing with load applied in the longitudinal 
direction, strain data were available from DIC and potentiometers. It was found 
that the measurements from potentiometers were more accurate, and therefore, 
this was used for evaluation of vlr and vlt- Poisson’s ratios for all tests are 
shown in Fig. 4. It can be seen that for certain directions, especially when 
evaluating very small strains in grain direction (due to high £), there is a large 
spread in results. In several of these cases (v TL and v RL ), measured displacements 
were less than 5 pixels, and although the software applies sub-pixel analysis 
techniques, this gave issues with accuracy. 




Fig. 2 Camera setup during experimental testing, a Compression test, b shear test 
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Fig. 3 Evaluation of Poisson’s ratio, a Definition of v, b grid on specimen, c horizontal and vertical 
strain and Poisson’s ratio 


Modulus of rigidity 

For calculation of the shear stress EN 408 (2010) specifies Eq. 3. This equation 
calculates the shear force and divides this over the full cross-sectional area of the 
specimen, assuming a constant shear distribution. This assumption has been checked 
using two-dimensional and three-dimensional FEM modelling. As there were only 
minor differences between 3D and 2D models, the 2D model (Fig. 5) is described here. 
For this model, steel support plates and timber have been modelled as generalised 
plain-strain elements. The epoxy connection between steel and timber has been 
modelled using tie constraints, a model with an adhesive layer included gave an almost 
identical shear stress distribution. Four node quadratic elements (CPEG8) have been 
used. Reduced integration elements (CPEG8R) were also tried and gave the same 
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Evaluation of the six Poisson’s ratios of LVL using digital image correlation software 
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results as fully integrated elements. One edge of the steel plate was loaded with a 
concentrated force of, respectively, 20 and 80 kN. This edge was restrained against 
movement perpendicular to direction of loading. The opposite steel plate had one 
restraint in direction of load and one perpendicular to direction of load. 



F max X COS 14° 

lb 



Figure 5 shows that shear stress distribution along the centre cross-section is not 
constant, but is highest in the middle 60 % of the specimen and drops to zero at both 
ends, which should be the case as there cannot be shear stresses at outside faces. 
From this model, it follows that maximum actual shear stress in the specimen is 1.22 
times the average shear stress. Therefore, when calculating the shear stress based on 
experimental testing this should be multiplied by a correction factor of 1.22. A very 
similar shear distribution has been found by Hassel et al. (2009). The almost uni¬ 
form shear stress distribution found by Sretenovic et al. (2004) could not be 
reproduced. 

The shear strain has been obtained using the DIC technique. Images were taken 
of each specimen at 5 second intervals until failure of the specimen. Using software, 
a grid of markers was placed over the centre section of the specimen (Fig. 6a) and 
markers were tracked through the series of images. An example of the displaced 
shape, multiplied with a factor 20, is shown in Fig. 6b. For each of the images, the 
horizontal and vertical shear deformations have been analysed. These deformations 
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Fig. 5 FEM model of shear test 
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Y = Yi + Y2 
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Fig. 6 Evaluation of shear strain using digital image correlation, a Specimen (Glt) with generated grid, 
b deformed grid (deformations multiplied by 20), c converting measured values into shear deformation 
and rigid body rotation 


have been added to get the total shear deformation (Fig. 6c), and rigid body rotation 
has been removed. When analysing this for all images the shear strain during testing 
can be found. Using timestamp of the images and timestamp of loading data, a 
stress-strain graph can be plotted. The slope of this graph, taken between 10 and 
40 % of maximum shear stress, gives the shear modulus. This procedure has been 
performed for all 36 shear tests, the result of which can be seen in Fig. 7. 
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Results and discussion 

The average and coefficient of variation values for the moduli of elasticity ( E ), 
Poisson’s ratios (v) and shear moduli (G) in the three material orientations are 
shown in Table 4. Values of E and G have been calculated using a lognormal 
distribution as outlined in AS/NZS 4357.3:2006 (Standards New Zealand 2006). 

As expected for timber, the stiffness values parallel to grain were much higher 
than the stiffness perpendicular to grain. The measured value of of 12.2 GPa is 
larger than the specified E of 11 GPa. This is to be expected as the E in compression 
is higher than the E in bending, and also because the E in bending includes shear 
deformation (Nelson Pine Industries Limited 2010). The stiffness in the tangential 
direction was slightly higher than in the radial direction. Variations in stiffness were 
less than 10 % for all three directions. It should be kept in mind that this is the 
material variability within one batch of LVL and is not representive of the 
variability within the LVL product. 

The Poisson’s ratio is strongly dependent on the direction, as is common for an 
orthotropic material. The variation in measurements is highest when extension 
happens in the longitudinal direction, vtl and vrl- This is due to small stresses in 
this direction and a high modulus of elasticity, resulting in very small strains. The 
second part of Table 4 shows average and coefficient of variation for off-diagonal 
terms in the compliance matrix. Normally, this matrix is assumed symmetrical, 
something which can now be verified. For that a t test has been used to determine if 
sets of two non-diagonal terms (independent samples with unequal variances) are 
statistically different. Although timber material properties are assumed to have a 
lognormal distribution, the skewness of the test data was sufficiently small to use a 
normal distribution. When comparing \’lt/El and v TL / E T the T value equals 0.16, 
for v lr /E l and v RL /E R the T value equals 0.79 and for v TR / E T and v RT / E R the T 
value equals 7.03. For most comparisons a total of 10 specimens were available, 
resulting in 8 degrees of freedom. The probabilities corresponding to the T values of 
0.16, 0.79 and 7.03 are 88, 45 and 0.01 %, respectively. This indicates a very good 
correspondence for the first test, a reasonable correspondence for the second test and 
a very certain difference for the last test. Bodig and Jayne (1982) have also verified 
symmetry of the compliance matrix. Their conclusion was that although deviations 
exist the assumption of symmetry holds reasonably well. Furthermore, they note 
that Vtl and v RR are nearly always quite small and may be subject to large 
experimental error, which has also been the case here. Although some difference 
between non-diagonal terms is measured, it would have large implications when 
implementing this in FEM software. Therefore, it is suggested to use Poisson’s 
ratios with the lowest variation: v R t, Vl R and vj R . 

The shear moduli found under six different loading orientations are listed in 
Table 4. It can be seen that there is a very good correspondence between Glt and 
Gtl- This should be the case as they correspond to loading in orthogonal plains 
which result in the same shear stress distribution and thus should result in the same 
shear stiffness. In addition, the rolling shear moduli, Gt R and G R t, are very similar. 
Only Gl R and G R l show a larger difference. In the final part of Table 4, 
experimental results of orthogonal plains have been averaged. These values show a 
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Fig. 7 Shear stress versus shear strain for six different directions 
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Table 4 Experimental and 
computed average values and 
coefficient of variations (CoV) 
for E, v and G 


Bold values are recommended 
for use in FEM software 


Property 

Average 

CoV (%) 

El 

12157 MPa 

7.3 

Ey 

426 MPa 

5.0 

Er 

371 MPa 

8.6 

'LT 

0.59 

32.4 

VTL 

0.02 

38.0 

'LR 

0.48 

17.6 

vrl 

0.02 

70.3 

Vtr 

0.22 

5.4 

VrT 

0.14 

12.8 

G\_t 

842 MPa 

14.3 

Gjl 

870 MPa 

11.8 

Glr 

1015 MPa 

8.3 

Grl 

786 MPa 

14.9 

Gjr 

103 MPa 

4.0 

Grt 

89 MPa 

10.8 

Vlt/El 

48 x 10 6 

29.4 

Vtl/ Ej 

47 x 10 6 

40.6 

Vlr / El 

39 x 10 6 

16.3 

Vrl! Er 

54 x 10 6 

70.9 

Vtr/Et 

516 x 10 6 

2.6 

Vrt/Er 

377 x 10 6 

10.8 

Glt+tl 

856 MPa 

13.2 

Glr+rl 

901 MPa 

18.0 

Gjr+rt 

96 MPa 

10.8 


variation between 10 and 18 %, indicating good repeatability in testing and analysis 
procedures. 

When comparing experimental results with literature values, shown in Table 1, 
several conclusions can be drawn. The modulus of elasticity in the longitudinal and 
the tangential direction (Er and Ej) are within range of literature data. Though in 
the radial direction (Er) the value is lower than all literature values. For sawn timber 
the modulus of elasticity in the radial direction is always larger than in tangential 
directions, but this does not seem to be the case for LVL. This was most likely due 
to stiffer and weaker layers sharing the load perpendicular to grain, whereas the 
weakest layers are governing the stiffness in the radial direction. On the contrast, for 
sawn timber stiffness in radial direction is higher than stiffness in tangential 
direction due to wood rays, which act as reinforcing rods in the radial direction 
(Wangaard 1979). Due to peeling and glueing of LVL, these rays are most likely 
damaged and do not contribute to the increased strength and stiffness in radial 
direction. The glue layers could also have contributed to the increased stiffness in 
tangential direction. Comparison of Er with other LVL products is not possible as 
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Janowiak et al. (2001) based modulus of elasticity values on five-point bending tests 
and thus could not evaluate the radial direction. 

The Poisson’s ratio vlt found by experimental testing of 0.59 corresponds well 
with values for other LVL products found by Janowiak et al. (2001), but within the 
literature data, there is a large spread in this Poisson’s ratio. The Poisson’s ratios vlt 
and v LR of Radiata Pine LVL are larger than of sawn pine species. The two small 
Poisson’s ratios, vtl and vrl, of 0.02 are within range of literature data. The values 
of vtr and especially vrt are below values mentioned in literature. This can be due 
to influence of gluelines which add stiffness and reduce strains in tangential 
direction. 

The shear moduli Glt and Glr of 856 and 901 MPa, respectively, are similar to 
values found in literature for softwoods. Shear stiffness values of different types of 
LVL, published by Janowiak et al. (2001), are lower than values found for Radiata 
Pine LVL and also mostly lower than other published values. For most species Glr 
is slightly higher than Glt, which is also found for Radiata Pine LVL. The rolling 
shear modulus, Gtr of 96 MPa is within range of values reported in literature. 


Conclusion 

An experimental testing campaign has provided elastic material properties for New 
Zealand Radiata Pine LVL. Block compression testing in three different material 
directions has resulted in moduli of elasticity of 12,157, 426 and 371 MPa in 
longitudinal, tangential and radial directions, respectively. It has been found that 
modulus of elasticity in the radial direction is lower than in the tangential direction, 
which is in contrast with literature values of sawn timber. The use of digital cameras 
and DIC has made it possible to extract Poisson’s ratios during the same tests. 
Poisson’s ratios compare well with values of other LVL products found in literature, 
but differ from literature values of sawn timber. Shear moduli have been determined 
using DIC on testing according to the European timber testing standard. Values of 
856 MPa and 901 MPa have been found for longitudinal shear moduli and 96 MPa 
for the rolling shear modulus. These results are in line with literature values for 
sawn timber, but higher than those for LVL manufactured from different species. 

The use of DSLR cameras for DIC has proven to give very satisfactory results, 
although accuracy for two Poisson’s ratios was limited due to very small strains 
along the grain direction whereby variations of up to 70 % were found. Variation in 
results for moduli of elasticity was less than 9 % and for shear moduli less than 
17 %. For most tests, only six specimens were tested, a larger sample size would 
improve the statistical basis. 

The reported elastic material properties make it possible to generate the 
constitutive matrix needed for FEM analysis. Diagonal terms follow directly from 
moduli of elasticity and shear moduli. Non-diagonal terms, calculated using 
Poisson’s ratios and moduli of elasticity, are not all symmetrical, but for practical 
purposes, assumptions can be made to keep the benefits of a symmetrical 
constitutive matrix. 
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